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Abstract 

We consider the motion of a spinning relativistic particle in ex- 
ternal electromagnetic and gravitational fields, to first order in the 
external field, but to an arbitrary order in spin. The noncovariant 
spin formalism is crucial for the correct description of the influence 
of the spin on the particle trajectory. We show that the true coordi- 
nate of a relativistic spinning particle is its naive, common coordinate 
r. Concrete calculations are performed up to second order in spin 
included. A simple derivation is presented for the gravitational spin- 
orbit and spin-spin interactions of a relativistic particle. We discuss 
the gravimagnetic moment (GM), a specific spin effect in general rela- 
tivity. It is shown that for the Kerr black hole the gravimagnetic ratio, 
i.e., the coefficient at the GM, equals unity (just as for the charged 
Kerr hole the gyromagnetic ratio equals two). The equations of mo- 
tion obtained for relativistic spinning particle in external gravitational 
field differ essentially from the Papapetrou equations. 
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1 Introduction 



The problem of the motion of a particle with internal angular momentum 
(spin) in an external field consists of two parts: the description of the spin 
precession and accounting for the spin influence on the trajectory of motion. 
To lowest nonvanishing order in c~ 2 the complete solution for the case of 
an external electromagnetic field was given more than 70 years ago [pj. The 
gyroscope precession in a centrally symmetric gravitational field had been 
considered to the same approximation even earlier |J . Then, much later the 
spin precession was investigated in the case of the gravitational spin-spin 
interaction M. The fully relativistic problem of the spin precession in an 
external electromagnetic field was also solved more than 70 years ago 
and then in a more convenient formalism, using the covariant vector of spin, 
in - 

The situation with the second part of the problem, which refers to how 
the spin influences the trajectory, is different. Covariant equations of motion 
for a relativistic spinning particle in an electromagnetic field were written 
in the same paper ||, and for the case of a gravitational field in ||. These 
equations have been discussed repeatedly from various points of view in nu- 
merous papers (see, e.g., [7-18]). The problem of the influence of the spin on 
the trajectory of a particle in external fields is not only of purely theoretical 
interest. It attracts attention being related to the description of the motion 



of ultrarelativistic particles in accelerators [[L9|] (see also recent review pOfl ). 

In fact, it is far from being obvious whether one can observe in prac- 
tice the discussed spin corrections to the equations of motion of elementary 
particles, say, electron or proton. According to the well-known argument 
by Bohr (see ||21||), an additional Lorentz force due to the finite size of the 
wave packet of a charged particle and to the uncertainty relation, exceeds 
the corresponding component of the Stern-Gerlach force. However, this ar- 
gument by itself does not exclude in principle the possibility to observe a 
regular Stern-Gerlach effect, let even a small one, in the presence of a larger 
background due to the uncertainty relation. This possibility is supported, in 
particular, by recent calculations | 22| . Moreover, spin-dependent correlations 
certainly exist in differential cross sections of scattering processes. So, it was 
proposed long ago to separate by polarizations a beam of charged particles in 
a storage ring through the spin interaction with external fields |[23|| . Though 
this proposal is being discussed rather actively (see review [p0|), it is not 
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clear up to now whether it is feasible technically. 

There are however macroscopic objects for which internal rotation cer- 
tainly influences their trajectories. We mean the motion of Kerr black holes 
in external gravitational fields. This problem is of importance in particu- 
lar for the calculation of the gravitational radiation of binary stars. In this 
connection it was considered in [24-27]. However, when we turned to these 
calculations, we found that the equations of motion with spin taken into 
account to the lowest nonvanishing order in c -2 , used in these papers, even 
in the simpler case of an external field lead to results which differ from the 
well-known gravitational spin-orbit interaction. The problem is essentially 
related to the correct definition of the center-of-mass coordinate. Moreover, 
it turned out that the widely used Papapetrou equations || also fail to repro- 
duce in the same c~ 2 approximation the result for the gravitational spin-orbit 
interaction found in the classical work . This discrepancy was pointed out 
long ago in f2"9fl ; however the explanation suggested in |29[ does not look 



satisfactory (see [p8fl ) 



In our recent work |3(J (the present talk is essentially based on it) we 
derived equations of motion of a relativistic particle with the noncovariant 
description of spin. These equations agree with well-known limiting cases. 
Though for external electromagnetic field such equations in the linear in spin 
approximation have been obtained previously |L9| (see also [0), we would 



like to start with comments related to this approximation in electrodynamics. 

2 Covariant and Noncovariant Equations 
of Motion of a Spinning Particle 
in an Electromagnetic Field 

2.1 The Problem with Covariant Equations of Motion 

The interaction of spin with external electromagnetic field is described, up 
to terms on the order of c -2 included, by the well-known Hamiltonian (see, 
for example, | |31|j ) 

H = -^sB+ e{9 - 1) s[pxE]. (1) 
2m 2m 2 w J V ; 
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Here B E are external magnetic and electric fields; e, m, s, and p are the 
particle charge, mass, spin, and momentum, respectively; g is its gyromag- 
netic ratio. Let us emphasize that the structure of the second (Thomas) term 
in this expression not only has been firmly established theoretically, but has 
also been confirmed with high accuracy experimentally, at any rate in atomic 
physics. To avoid misunderstandings, let us note that, generally speaking, 
the last term in formula ([!]) should be rewritten in a Hermitian form (see, 
e.g., II): 

[px E}^^ {[px E\-[E xpj) = [px E]+ l -V x E. 

We will be interested, however, in the semiclassical approximation mainly, 
when in the interaction linear in spin, field derivatives are neglected. (Be- 
sides, the correction with V x E vanishes in the case of potential electric field 
considered in [|3"I[].) 

Let us try to construct a covariant equation of motion accounting for spin, 
which would reproduce in the same approximation the force 

r - = ^ gS - + (jt [B x i- " x e - ] ) • (2) 

corresponding to the Hamiltonian ([I]) (here and below coma with a subscript 
denotes a partial derivative). A covariant correction to the Lorentz force 
eF^ v u vl linear in the tensor of spin and in the gradient of the tensor 
of electromagnetic field F^x , may depend also on the 4-velocity Since 
u^u^ = 1, this correction must satisfy the condition u^f 1 = 0. From the 
mentioned tensors one can construct only two independent structures meet- 
ing the last condition. The first, 

rTFuX^ - F Xu , K u*S x »u^ (3) 

reduces in the c -2 approximation to 

2s(B, m -[v x E, m ]), 

and the second, 

u x F Xv!K u K S^, (4) 
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reduces to 

(Let us note that the structures with the contraction F vk \S kX reduce to these 
two expressions, due to the Maxwell equations and the antisymmetry of S K \.) 

Obviously, no linear combination of these two structures can reproduce 
the correct expression (0) for the spin-dependent force. In a somewhat less 



general way it was shown in [28 



But why the correct (in the c -2 approximation) formula (|2|) cannot be 
obtained from a covariant expression for the force? Obviously, one can easily 
reproduce those terms in (0) which are proportional to g by a linear combina- 
tion of (§) and (U), i.e., there is no problem to present in a covariant form the 
terms which describe, so to say, direct interaction of a magnetic moment with 
external fields. It is the terms in (|2|) independent of g and corresponding to 
the Thomas precession which cannot be written covariantly. Certainly, the 
Thomas precession can be described beyond the c~ 2 approximation, for ar- 
bitrary velocities. But there are no reasons why this essentially noncovariant 
phenomenon should have a covariant description. This is the point. 

2.2 What Is the Correct Definition of the Coordinate 
of a Spinning Particle? 



It was noted in [g8| that the covariant formalism can be reconciled with the 
correct results if the coordinate x entering the covariant equation is related 
to the usual one f in the c -2 approximation as follows: 

1 _ _ 

x = r + — s x v. (5) 
2m 

. The generalization of this substitution to the case of arbitrary velocities 

7 1 
x = f+ — — — s xv, 7 = (6) 

m(j + l) \/l-v 2 



was recently pointed out in p0 | 



But why the spin precession itself (as distinct from the spin influence on 
the trajectory) can be described covarianly [||, [5] without any concern for 



4 



the coordinate definition? First of all, the covariant equations of the spin 
precession 

dSfj, e 



(h v gF^ -(g- 2) Ufl F Xu u x S»] (7) 

(here is the covariant 4-vector of spin) are written in the semiclassical 
approximation, i.e., the coordinate dependence of external fields is neglected 
at all. Second, equations (|7|) are homogeneous and linear in spin. So, even 
if one went here beyond the semiclassical approach, but stayed within the 
approximation linear in spin, the use of the usual coordinate r, which differs 
from x in terms proportional to s only, would be completely legitimate. 

Of course, the choice of the variable, r or x, is by itself a matter of 
convention. But still, which of them is the true center-of-mass coordinate of 
a relativistic spinning body? 

We note first of all that relation (^) is valid for a free particle as well. So, 
to answer the question, we will consider a simple example of the free Dirac 
particle with the Hamiltonian 

H = ap + j3m . 

Here, the operator whose expectation value equals to r, is not r itself, but |33| 

2 = f+ <g* _ ^N, ; e = g = i [Sxa1 . (8) 

To lowest nonvanishing order in c~ 2 expression (^) reduces to 

1 ^ - - 1^ 
x = r + s x v , s = -a , (9) 

which might prompt indeed substitution ([5]). However, under the Foldy- 
Wouthuysen (FW) transformation, which separates positive-energy states 
from negative-energy ones, the relativistic operator x goes over into mere 
r. And the transition from the exact Dirac equation in an external field 
to its approximate form containing only the first-order correction in c~ 2 is 
performed just by means of the FW transformation. Thus, in the arising 
c~ 2 Hamiltonian the coordinate of a spinning electron is the same r as in 
the completely nonrelativistic case. Nobody makes substitution ([5]) in the 
Coulomb potential when treating the spin-orbit interaction in the hydrogen 
atom. 
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As to a classical particle, it is in fact a well-localized wave packet con- 
structed from positive-energy states, i.e., it is properly described in the FW 
representation. Therefore, it is r which is the true coordinate of a classical 
relativistic spinning particle. 



2.3 The Noncovariant Formalism 



The correct equations of motion in an electromagnetic field including spin 
to first order are known for a fairly long time Though being fully 

relativistic, they are noncovariant and based on the initial physical definition 
of spin. According to this definition, spin is the 3-dimensional vector s (or 
3-dimensional antisymmetric tensor s mn ) of the internal angular momentum 
defined in the rest frame of the particle. The covariant vector of spin S ^ 
(or the covariant antisymmetric tensor S^) are obtained from s (or s mn ) 
merely by the Lorentz transformation. By the way, an advantage of this 
approach is that the constraints u^S^ = and u^S^ = hold identically. 
The precession frequency for spin s at an arbitrary velocity is well-known 
(see, for instance, [BID: 



B 



7 



7 + 1 



v(vB) 



v x E 



1 - 1 -> 1 

+2 -B v x E } . (10) 

Naturally, the corresponding interaction Lagrangian (the Lagrangian descrip- 
tion is here somewhat more convenient than the Hamiltonian one) equals 



L u = Qs= —s Ug - 2) 



B 



7 



7 + 1 



v(vB) 



v x E 



+2 



-B — v x E 

7 7 + 1 



(11) 



The equation of motion for coordinate is the usual one: 



(12) 
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where L to t is the total Lagrangian of the system. The equation of motion for 
spin in general form is 

s = -{L tot ,s}, (13) 
where {...,...} is the Poisson bracket, or 

k=-i[L tot ,s\ (14) 

in the quantum problem. This applicability of a common canonical formalism 
is one more advantage of the noncovariant approach. Meanwhile, in the 
covariant approach we have to deal in particular with higher time derivatives, 
which is obvious already from relation (|J). 



3 Equations of Motion of a Spinning Particle 
in an Electromagnetic Field 

3.1 The General Approach 

Our approach to the problem is based on the following physically obvious 
argument. As long as we do not consider excitations of internal degrees 
of freedom of a body moving in an external field, this body (even if it is 
macroscopic) can be treated as an elementary particle with spin. 

Therefore, the Lagrangian of the spin interaction with an external field 
can be derived from the elastic scattering amplitude 

- eJ^, (15) 

of a particle with spin s by a vector potential A^. Due to the arguments pre- 
sented in Section 1, the discussion of the effects nonlinear in spin (which are 
of primary interest to us) may be physically meaningful only in the classical 
limit s ^> 1. It is this approximation that is basically used below. 

The matrix element of the electromagnetic current operator between 
states with momenta k and k! can be written (under P and T invariance) as 



follows (see, for instance, 34, K 



Here p M = (k' + k)^ g M = (k' - k)^ 
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The wave function of a particle with an arbitrary spin ijj can be written 




(see, for instance, ]3T|, §31) as 

^ = -$s ( ; ) • (17) 

Both spinors, 
and 



$ ~ ^$102- Pq 



7 I'/ai a?2 •■ a P J ' 

are symmetric in the dotted and undotted indices separately, and 

p + q = 2s. 

For a particle of half-integer spin one can choose 

1 1 

P = s+ -, q = s- - . 

In the case of integer spin it is convenient to use 

p = q = s. 

The spinors £ and rj are chosen in such a way that under reflection they 
transform into each other (to within a phase). For p ^ q they are differ- 
ent objects which belong to different representations of the Lorentz group. 
For p = q, these two spinors coincide. Nevertheless, we will use the same 
expression ([T7|) for the wave function of any spin, i.e., we will also formally 
introduce the object rj for an integer spin, bearing in mind that it is ex- 
pressed in terms of £ . This will allow us to perform calculations in the same 
way for the integer and half-integer spins. 

In the rest frame both £ and 77 coincide with a nonrelativistic spinor £0, 
which is symmetric in all indices; in this frame there is no difference between 
dotted and undotted indices. The spinors £ and 77 are obtained from £ 
through the Lorentz transformation: 

f = exp{S0/2K o ; V = exp{-S0/2}£ o • (18) 
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Here the vector is directed along the velocity, tanh0 = v\ 

^ v p+q 

i=l i=p+l 

and &i acts on the ith index of the spinor £ as follows: 

^i& = (*i)a 4 ft (&)....&... ■ (19) 



In the Lorentz transformation ( IB ) for £, after the operator S has acted on £ 
the first p indices are identified with the upper undotted indices and the next 
q indices are identified with the lower dotted indices. The inverse situation 
takes place for r\. 

We note that in an external field the components of velocity v (and to- 
gether with them the components of 0) do not commute, in general. However, 
to the adopted approximation, linear in the external field, one can ignore this 
noncommutativity which is itself proportional to the field. Moreover, we are 
mainly interested in the classical limit of the final result where such com- 
mutators are negligible since they are proportional to an extra power of h. 
Therefore, v and <ft will be treated as ordinary numerical parameters. 

Next, 

o l\ 

/or 



■0 = V^To = 4> ] 



here / is the sum of 2 x 2 unit matrices acting on all indices of the spinors 
£ and rj . The components of the matrix = — are: 



J 0n 




(20) 

Emn = 2ze m?1 fc [ n ] i (^-Q 



The scalar operators F e ^ m depend on two invariants, t = q 2 and r = 
(S^q^) 2 . The covariant vector of spin is defined, e.g., for the state with 
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momentum k^, and is obtained via the Lorentz transformation from the 
vector of spin (0, s) in the rest frame: 

S"={S ,§), S = { -^-, S = s+ (22) 

m m(e + m) 

In the expansion in the electric multipoles 

N e 

F e (t,T) = Y J fe,2n(t)T n 

71=0 

the highest power N e is obviously s and s — 1/2 for integer and half- integer 
spin, respectively. In the magnetic multipole expansion 

N m 

F m (t,T) = J2fm,2n(t)r n 
n=0 

the highest power N m is s — 1 and s — 1/2 for integer and half- integer spin. 
The invariant form factors f e , m (t), being Fourier-transformed into the coor- 
dinate representation, describe the space distribution of the charge, magnetic 
moment, and higher multipoles, i.e., describe the finite size of the particle. 
Clearly, 

/ e ,o(0) = 1, / m , (0) = |. 

Let us note at last that we have chosen the noncovariant normalization 
for the amplitude (|T5|), being interested in the Lagrangian referring to the 
world time t and not to the proper time r. 

3.2 Effects Linear in the Spin 

As a warm-up exercise and check of our approach, let us reproduce now the 
well-known result (|TT]) for the case of a constant external field. We start with 
the terms proportional to ^-factor. The corresponding term in the scattering 
amplitude can be written as 

^ £o f {[exp{S0/2}(s5) exp{-S0/2} + exp{-E0/2}(£B) exp{S0/2}] 

+ l - [exp{S0/2}(Si) exp{-S0/2} (23) 
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- exp{-£0/2}(££) exp{£0/2}] } £ . 

It is essential that in the considered case of a constant external field, one 
may put k=k, v=v, (/)=(/), since q = k — k corresponds to the field 
gradient. 

In our further calculations we use the well-known identity 



exp{A}Bexp{~A} = B + 1 [i, 5] + 1 [i, [i, 5] 



+ 



and the following relationships: 

cosh = 7, sinh = t> 7. 
After simple algebraic transformations, expression (|23|) reduces to 



eg 
2m 



5- 



7 



7 + 1 



v(vB) — v x E 



Let us now discuss the contribution of the convection term 
We write the product of exponents in the expression 

= -$[exp{S^/2}exp{-E?/2} + exp{-£0/2} exp{£0/2}]£ o 



as 



(24) 
(25) 

(26) 
(27) 



(28) 



exp{E0 /2}exp{-E0/2} 
= ]Jexp{a<p /2}exp{-a<p/2}]]_exp{-a<p /2}exp{a<p/2}. (29) 

Let us consider a typical factor in this formula: 

exp{a0 /2} exp{-o$/2} = cosh(0'/2) cosh(0/2) - (nn) sinh(0'/2) sinh(0/2) 



+ a 



rasinh(072) cosh(0/2) - ncosh(0'/2) sinh(0/2) 



(30) 
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- i (a[n x n\) sinh(0'/2) sinh(0/2); 

here n = v /v', n = v/v . Since we are interested in gradients only as long 
as they enter together with spin, in the first term, cosh(0'/2) cosh(0/2) — 
(nn) sinh(0'/2) sinh(0/2), we put <$' = <f)/2, n = n , after which this term 
turns to unity. Then, we are discussing the interaction linear in spin, so that 
the product (|2~9D reduces to 



1 + £ [nsinh(072) cosh(0/2) - ncosh(0'/2) sinh(0/2) 
- 2i (s[n x n]) sinh(0'/2) sinh(0/2). 



When substituted into formula (0), the terms proportional to £ cancel out. 
Now, limiting ourselves to the terms linear in q, we reduce the spin-dependent 
part of (|7|) to 

p^ i(s[k x q\) . 
2e m(e + m) 

Let us note further that since p^q^ = , the following identity holds 

p^q a A^ = p^iqaA^ - q^A a ) = pHF^. (31) 

Then, we can put p^ = 2mu il , where is the 4-velocity. As a result we 
arrive at the following expression: 



2m 



2 1- 



B 



7, 



27 
7 + 1 



v(vB) 



27 
7 + 1 



v x E 



(32) 



The sum of fl26|) and ( |32"D yields (p+|). Thus, we have reproduced the well- 
known result for the interaction linear in the spin, starting from the relativis- 
tic wave equation for an arbitrary spin. 

Below we repeatedly use identities of the form ([31]). In the classical 
language such a transformation corresponds to discarding in a Lagrangian 
(or adding to it) a total time derivative. Indeed, 



7 



' d_ 
dt 



+ vV 



7 



d_ 
dt' 
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3.3 Effects Quadratic in the Spin 



Let us now investigate the interaction of second order in the spin. The "bare" , 
explicitly quadrupole interaction present in expressions (|i~5|) and (|16|) is 



(33) 



Using the identity ( [31] ) and relations (^), and then discarding the total time 
derivative d/dt + vV, we write this interaction as 



J 2s 



ef, 



e,2 



sE) ^_{sv){vE) + (s[vx B}) 



7 



7 + 1 



(vs){vV) 



7 + 



(34) 



Using the Maxwell equations and adding a total derivative with respect to 
t, one can show that the tensor s^Sj in (|34"D can be rewritten in the following 

SiSj — (l/3)5ijS 2 . Now from the nonrelativistic limit 



irreducible form: s^Sj 



of formula (p4[) , it is clear that this formula describes indeed the interaction 
with an external field of the quadrupole moment 

Qij = -2e fe^iSsiSj - 5ijS 2 ); Q = Q zz \ Sz=s = - 2 e f Ct2 s(2s - 1). (35) 

In the asymptotics, as 7 — > 00, the interaction (^) tends to a constant 

L 2s = -e/ e , 2 [(sV) - (vs)(vV)} [(sE) - (sv)(vE) + (s[v x B])]. (36) 

It is well-known that even in the absence of the bare quadrupole term, 
i.e., at f e ^2 = 0, a quadrupole interaction arises in the nonrelativistic limit 
due to the convection and magnetic terms in interaction (|T5|). The value of 
this induced quadrupole moment at an arbitrary spin of the particle is [|35| 
(in the formula below we corrected a misprint in the original paper ||35|| ) : 



Q 1= -e(g-l] 



fi \ j s, integer spin, 

mc J \ s — 1/2, half-integer spin. 



(37) 



Here we have explicitly displayed the Planck constant h to show that the 
induced quadrupole moment Qi vanishes in the classical limit h —>■ , 
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s — > oo, hs — » const. Therefore, the interaction of second order in spin 
proportional to Q\ does not influence in fact equations of motion of a classical 
particle (although it plays a role in atomic spectroscopy P5| ). 

Still, the convection and magnetic terms in expression (JT5]) induce an 
interaction of second order in spin which has a classical limit and is therefore 
of interest for our problem. It is convenient here to start with the convection 
current interaction. Let us come back to formula (130). Again we put in it 



cosh(07 2 ) cosh(0/2) - (nn) sinh(0'/2) sinh(0/2) = 1. 

In the other terms, linear in a, we keep only the first power of q — > — iftV, in 
the hope that in the final result ( |29| ) % will enter in the combination hs — > 
const. Nevertheless, these terms by themselves are small as compared to 
unity, so that in the classical limit expression ( j30|) can be rewritten as 

exp jcr n'sinh(0'/2) cosh(0/2) — ncosh(0'/2) sinh(0/2) 

-i(a[n' x n])sinh 2 (0/2)}. 
Clearly, in the product the operators a attached to 

n sinh(0'/2) cosh(0/2) - ncosh(0'/2) sinh(0/2) , 

combine in the resulting exponent into the operator S which vanishes in the 
classical limit. In this limit only those operators a survive that are attached 
to [n x n] sinh 2 (0/2); they combine into 2s. Thus, in the classical limit the 
product ( p9|) reduces, with the account for the second identity ( p5|) , to 

expjl-^^xV])}. (38) 

Let us note that the action of the operator ( |38| ) on any function of co- 
ordinates, whether it is a vector potential or field strength, amounts to the 
shift of its argument: 

1 7 - 

r — > r H s x v . 

m 7 + 1 



Curiously, just this substitution was pointed out in [5D| for the transition 



from covariant equations linear in spin to noncovariant equations (see (0)). 
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Now, taking into account the second term in the expansion of the ex- 
ponential function fl38|) and using again the identity (0), we obtain the 
following expression for the quadratic in spin interaction arising from the 
convection current: 



7 



2m 2 7 + 1 



[s[vx V]) 



i - i m 



7 



(sv)(vB) — (s[vx E] 

7 + 1 v 



(39) 

7 + 1 7 + 

Let us discuss now the contribution to the discussed effect due to the mag- 
netic moment. It is convenient to write the term in the scattering amplitude 
we are interested in (it is proportional to g-factor) as 

^ £o f |[exp{E0/2}exp{-E0/2}exp{E0/2}(sB)exp{-E0/2} 



+ exp{-£0 /2} exp{£0/2} exp{-E0/2}(£B) exp{E0/2}] 
+ ~ [exp{E0/2} exp{-E0/2} exp{E0/2}(EE) exp{-E0/2} 

- exp{-S0 /2} exp{E<?/2} exp{-E0/2}(EE) exp{E0/2}]l $ 



(40) 



Using in this case the first term in the expansion of the exponential function 
(p8[) , we arrive at the following expression for the contribution proportional 
to the magnetic moment: 



(s[vx V]) 
2m 2 7 + 1 v L iJ 



(SB) 



7 



7 + 



-(3S7)(tfB) - (s[vx E) 



(41) 



The total result for the induced interaction, quadratic in the spin, is 



7 



2m? 7 + 1 



[s[v x V]) 



g - 1 + - ) (sB) -(g-1) -Jt (S?0(tfB) 
v 7/ 7 + 1 



7 



7 + 



(42) 



Let us note that in the nonrelativistic limit the induced interaction with 
magnetic field tends to zero as v/c, and that with electric field as (v/c) 2 . 
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Moreover, the part of interaction ( |4*2| ) that is not related to g-factor, is re- 



ducible in spin; in other words, SiSj in it cannot be rewritten as an irre- 
ducible tensor SiSj — (l/3)5ijS 2 . Therefore, the interaction ( f4"2| ) is not in 



fact a quadrupole one. However, its asymptotic behaviour for 7 — > 00 is of 
interest. In this limit 

LL = ^ (9 ~ 1) (s[v x V]) [(sB) - (sv)(vB) - (s[v x E})}. (43) 

Surprisingly, the asymptotical formulae (|36|) and ([|3|) coincide to within a 
factor and a total time derivative. To prove this, it is convenient to introduce 
three orthogonal unit vectors 

[vxs\ T 
\[v X s\\ 

Using the completeness of this basis and the equation E = [V X B] , and 
discarding a total derivative with respect to t , one can check that 

[(sV) - (vs)(vV)\ [(sE) - (mJ)(vE) + (s[v x B})} 
= [v x S f (CV)[(C^) + (P-B)], 

coincides indeed with 

x V]) [(sS) - (sv)(vB) - (s[vxE])\ 

= -[vxs\ 2 (pV) [(^ x S]) + (pE) . 
Thus, there is a special value of the bare quadrupole moment 

es 2 1 
Q = -2(g-l) — , or / e , 2 = (9-1)^ (44) 

(let us recall that we consider now a classical situation, when s ^> 1 ), at 
which the total interaction quadratic in spin, L 2s + L % 2s , asymptotically de- 
creases with energy. 

The situation resembles that which takes place for the interaction linear 
in spin. It is well-known (see, for instance, |TT|, |36|, 0) that there is a special 
value of ^-factor, g — 2 , at which the interaction linear in spin decreases 
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as 7 — > oo . This follows immediately from formula (|11|) for the first-order 
Lagrangian. Thus, putting additionally g = 2 , we obtain 



Q 




2m 2 



1 



(45) 



Let us note that the choice g = 2 for the bare magnetic moment is a 
necessary (but insufficient) condition of renormalizability in quantum elec- 
trodynamics [Tl], |36], [37]] . It is satisfied not only for the electron, but also for 
the charged vector boson in the renormalizable electroweak theory. 

In some respect, however, the situation with the special values (^4|), (fl5|) 
of the quadrupole moment differs from the situation with ^-factor. The 
conditions fl44]) , (fi5|), as distinct from the condition g = 2 , are not universal, 
since they are valid only for large spins, s ^> 1 ; in other words, they refer 
only to classical objects with internal angular momentum. In particular, for 
the charged vector boson of the renormalizable electroweak theory the bare 
quadrupole interaction is absent at all, f e ^ = . The quadrupole moment of 
this particle is (in our language) of the induced nature, it is given by formula 



In this section we present a simple and general derivation of the equations 
of spin precession in a gravitational field. This approach not only allows 
us to easily reproduce and generalize known results for spin effects. Pointed 
out here remarkable analogy between gravitational and electromagnetic fields 
allows also to easily transform the results of the previous section to the case 
of an external gravitational field. 

It follows from the angular momentum conservation in flat space-time 
taken together with the equivalence principle that the 4-vector of spin S 1 * is 
parallel transported along the particle world-line. The parallel transport of 
a vector along a geodesic x M (r) means that its covariant derivative vanishes: 



(|7P at s = 1 g = 2 . 



4 A Simple- Minded Aside on the Spin 
Precession in a Gravitational Field 



DS" 



(46) 



0. 



Dt 
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(In this section we restrict our discussion to the effects linear in the spin.) 
We will use the tetrad formalism natural for the description of spin. In view 
of relation (P>|), the equation for the tetrad components of spin S a = S^e^ 



is 

^=^f= S KX = V ab J bcd u d S c . (47) 

Here 

7abc = e atl . v e%e v c = -% ac (48) 



are the Ricci rotation coefficients Certainly, the equation for the tetrad 
4- velocity components is exactly the same: 

dn a 

• (49) 



dr 



The meaning of Eqs. (|4"TD, fl49| ) is clear: the tetrad components of both vectors 
vary in the same way, due to the rotation of the local Lorentz vierbein only. 

In exactly the same way, the 4-dimensional spin and velocity of a charged 
particle with the gyromagnetic ratio g = 2 precess with the same angular 
velocity in an external electromagnetic field, by virtue of equation (|?p at 
g = 2 and the Lorentz equation: 

dSg e b _ dua e b 

—, — — — PabO i —T~ — — fabU ■ 

dr m dr m 

Thus, the correspondence: 

g 

— Fab * ► labcU C . (50) 

m 

gets obvious. This correspondence allows one to obtain the precession fre- 
quency u of the 3-dimensional vector of spin s in external gravitational field 
from expression (|ToD via the simple substitution 



e l 6 

— Bi — > - -e ik ilk\cU c \ —Ei — > ^oi C u c . (51) 
ml m 



This frequency is 



Ui = -e iH -7h c H 7— —rime -fr • (52) 

V 2 u v + 1 / K„ 
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The factor l/vP w in this expression is related to the transition in the left-hand 
side of Eq. (|47[) to the differentiation with respect to the world time t: 



d dt d q d 
dr dr dt w dt 

The quantity w° is supplied with the subscript w to emphasize that this 
is a world component of 4- velocity, but not a tetrad component. All other 
indices in expression ( p)2"D are tetrad ones, c = 0, 1, 2, 3; i,k,l = 1, 2, 3. The 
corresponding spin-dependent correction to the Lagrangian is 

L lsg = suj . (53) 

As an illustration of formulae (0), Q53"|), let us apply them to the cases 
of spin-orbit and spin-spin interactions. We restrict, as is common in the 
problems discussed, to the linear approximation in the gravitational field. 
However, in our approach, as distinct from the standard ones, both problems 
can be easily solved for arbitrary particle velocities. 

The tetrads e aAt are related to the metric as follows: 

To linear approximation we can put g^ u = r] nv + h nv and do not distinguish 
between the tetrad and world indices in e afl . The ambiguity in the choice of 
tetrads will be fixed by choosing the symmetric gauge e nu = e vn . Then 



Using expression fl48 ) for the Ricci coefficients, we find to linear approxima- 
tion | 

labc = ^(hc,a- Kc,b) ■ (54) 

Let us start with the spin-orbit interaction. In the centrally symmetric 
field created by a mass M, the metric is 

2kM 2kM 

^00 — — i ""mn — _ — mn . {00) 

Here the nonvanishing Ricci coefficients are 

lijk = -^-{OjkTi - 6 ik rj) , 7oi = — ^T r i- I 56 ) 
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Their substitution into formula ( [52] ) yields the following expression for the 
precession frequency: 



2 7 +lfcM^ _ 

uis = — — v x r . (57) 

7+1 r- 5 

In the limit of low velocities, 7 — > 1, the answer goes over into the classical 
result @. 

Now we consider the spin-spin interaction. Let the spin of the central 
body be sq. Linear in sq components of metric, which are responsible for the 
spin-spin interaction, are: 

h x = 2k [? ° X 1' 



Here the nonvanishing Ricci coefficients are 

7*=*(vK^-V,K£3i), 7W ^-,V,I^. (58) 
The frequency of the spin-spin precession is 

uj.. = —k 12- 




+ k [v(s V) - s (vV) + (vs )V] (tfV) - . (59) 

7+1 r 

In the limit of low velocities this formula also goes over into the corresponding 
classical result 0. 

In the conclusion of this section we note that in the case of an external 
gravitational field there is no covariant expression for the force linear in the 
particle spin. In other words, the deviation from geodesies of the trajectory 
of a spinning particle is not described by the Riemann tensor. If it were 
the case, there would be a unique possible covariant structure, to within 
a factor (in it equals — l/2m): R^ uab u v S ab . As mentioned already in 
Section 1, the covariant description (as distinct from our our formulae ([52]) , 
contradicts the classical results in the limit of low velocities. 
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5 Equations of Motion of a Spinning Particle 
in a Gravitational Field 

5.1 The General Approach 

The equations of motion in an external gravitational field to any order in 
spin are constructed similarly to the equations of motion in the case of an 
electromagnetic field. 

We start with the elastic scattering amplitude in a weak external grav- 
itational field hpv. We use it as an euristic argument only, and afterwards 
will go beyond the linear approximation. This amplitude is 

- \T^ V . (60) 

The matrix element T^ v of the energy-momentum tensor between states of 
momenta k and k! can be written as 

= — 4){k') Ip^pv Fi + i (p^ uX + Pu^^x) q X F 2 

+ (v^q 2 ~ %qv) F 3 (61) 

+ [S^q 2 - (5^ + S vq[i )(Sq) + V^(Sq) 2 ]F A } ^(k). 

The scalar operators in this expression are also expanded in powers of 
r = (Sqf: 

Ni 

F i (t,r) = X;/i,2n(t)r n . 

n=0 

It can be easily proven that the total number of invariant form factors f^n 
is 4s + 2 and 4s + 1 for integer and half-integer spin, respectively. The 
independence of the four tensor structures in ( |6T]) is obvious. As for the 
completeness of the expansion, it can be proven, for instance, by demon- 
strating that the total number of invariant form factors, as calculated in the 
annihilation channel, coincides with the above result. 

In the generally covariant form, the structure (rj^q 2 — q^q v ) W cor- 
responds to the scalar curvature R, and [S^S v q 2 — (S^q v + S v q^)(Sq) + 
Tj^ v {Sq) 2 } h^ u corresponds to the product R^S^S" ', where R^ u is the Ricci 
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tensor. Since we are interested in the equations of motion in a sourceless 
field, the corresponding terms in the expansion ([[]]) will be omitted. 

Just as in electrodynamics charge conservation dictates the condition 
/e,o(0) = 1 , here energy conservation leads to /i,o(0) = 1 • As to the term 
in the amplitude (|60f) which contains /2,o > h is convenient to write it in a 
different form, using the analogy fl50|) with electromagnetic field. Putting 
g = 2 and (e/m)F ab — > f ab = -y a b c u c in the corresponding electromagnetic 
term 

i^-$(k') Z ab F ab ^(k), 
oe 

we arrive at the following contribution to the Lagrangian of the gravitational 
interaction: 

i^^(k')E ab f ab ^(k); (62) 

here, as usual, u° = e/m. Using for •y a b c the linear approximation (|54]), one 
can easily see that the expression (|62"D corresponds indeed to the discussed 
contribution to the amplitude provided that / 2 ,o — 1- Thus, in gravity 
the value of one more form factor at zero momentum transfer t is fixed. It 
corresponds to conservation of angular momentum. This was pointed out 



earlier in fl39l 40 



Let us now come back to the convection term in formula (0). As in 
electrodynamics, when we go over here to spinors in the rest frame, the term 
of first order in spin is written as 

- ^-irr(^ x V]) K v . (63) 
Using (|31~D, (p4|), we obtain 

P^kh^u -> -P^i-dkh^ + d^h kv ) -> -2p a >y aku . 
Thus, the expression (|63|) can be written in terms of the Ricci coefficients: 

1 7/° 

- — ^Wr*^. (64) 



As can be easily seen, the sum of (|^) and ( |64D reproduces the Lagrangian 
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5.2 Effects Quadratic in the Spin 

Let us investigate now the effects of second order in the spin in the equations 
of motion in a gravitational field. In the case of a binary star these effects 
are of the same order of magnitude as the spin-spin interaction when the 
spins of the components of the system are comparable [p8 |. The influence of 
the latter on the characteristics of gravitational radiation becomes noticeable 
for a system of two extreme black holes [p5fl . Correspondingly, second-order 
spin effects in the equations of motion become substantial if at least one 
component of a binary is close to an extreme black hole | 28| . Therefore, the 
investigation of these effects is not of a purely theoretical interest only. In 
principle they can be observed with the gravitational wave detectors under 
construction. 

An obvious source of second-order spin effects is the term 

L2s 9 = -A,2^PW<?)V (65) 
in the amplitude fl60|). Due to the relation 

p^p u q a qph^ = p fx p u {q a q f 3h^ + q^Kp - q a q u h^ - qpq^h ua ) 

-> 2^p v R^, 

the Lagrangian (|65|) can be written in terms of the Riemann tensor: 

i2w= _^ S VrtW (66) 
Instead of /i j2 , we have introduced a dimensionless parameter k: 

/l2= — . 

2m 2 

Now, it is convenient to use the Petrov representation for the components 
of the Riemann tesor (see 



Ekl — RokOl, Ekl — Elk] Ckl — -^kmn^-lrsRrnnrs j CfcZ — C[k] 

B U = -e lrs R 0krs , B kk = 0. (67) 
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We limit our discussion to the case of a sourceless gravitational field. Then, 
at R a b = 0, further simplifications hold: 



C, 



E, 



kh 



Bki = Bik, E, 



ik- 



J kk 



Ckk — o. 



(68) 



Finally, we arrive at the following interaction Lagrangian quadratic in the 
spin: 



K 

Ye 



(2u 2 + 1) E kl - 2 
1 

+ 



(u + iy 



u + 1 

u k uiu m u n E, t 



UkU m Elm + 5klU m U n E 

II 



(69) 



2 uq e}~ mri u m B n i 4 



V>kU m €i rn U r B mn 



[s k si - - 5 M s 2 ) 



u + 1 

To avoid misunderstandings, we note that all three-dimensional indices in 
this equation (and in ([70|), (|7l|) ) are in fact contravariant. 

As in electrodynamics, along with the "bare" interaction (69), there is 
here an induced interaction quadratic in spin. Its explicit form can be ob- 
tained most easily by setting g = 2 in the electromagnetic formula (E2J) and 
by making the substitution (|5lD. We also take into account the correspon- 
dence 



qilabcU = (qaabc ~ qdabi) u c -> % {dcj ahc - d c ^ ab i) u c 



^ ELabci'U • 



Finally, using ( p7|) and (|6q) , we obtain the following result for the induced 
interaction: 



J 2sg 



1 
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2vf 



u 



u°+ 1 

1 



E, 



ki 



1 



1 



-2 k 



(«o + I) 2 
1 

U + 1 



8klU m UnE n 
^■kmn'U'mBnl ~\" 



+ 



U + 1 
1 

(uo + 1) : 



(«o + I) 2 

u k uiu m u n E, 



UkUmElm 
,n (TO) 



UkUrn^lrnUrB mn f S^Sl- 



Uq + 1 

As in the electromagnetic case, the induced interaction tends to zero in the 
nonrelativistic limit ~ v/c, and the spin factor in it, SkSi , is not an irre- 
ducible tensor. 
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The asymptotic behaviour of L2 Sg and L\ sg is the same: both Lagrangians 
increase linearly with energy. However, in this case too the coefficient in 
the "bare" interaction can be chosen in such a way, k — 1 , that the total 
Lagrangian of second order in the spin decreases (as well as the analogous 
interaction in electrodynamics) when the energy tends to infinity. At k = 1 

L 2sg + Li S g = -^o^i) (u°E kl - -^- iUk u m E lm 

SkSi- a 

2(u +1) J 

6 The Gravimagnetic Moment 

There is a profound analogy between the linear in the spin interaction of the 
magnetic moment with the electromagnetic field 

£ e m = —-r-F a i 1 S ab (72) 
4m 



and the "bare" gravitational Lagrangian (|6q) which is quadratic in the spin [|TT 
(Here it is more convenient to write the gravitational Lagrangian, like C em , 
for the proper time r , rather than the world time t , i.e., to multiply expres- 
sion (|66|) by e/m.) This analogy is based on the following observation [[11 



It is well-known that the canonical momentum id^ enters relativistic wave 
equations for a particle in electromagnetic and gravitational external fields 
via the combination 

n M = zfy - eA^ - h: ab iab». 

It follows from the structure of the commutator (or Poisson brackets in the 
classical limit) 

1. 



[n^n,] = -i(eF^- -Y> ab R ab ^) 

that in a sense — |S afc i? abAti/ plays the same role in gravity as eF^ u in elec- 
tromagnetism. It is quite natural then that the gravitational analogue of the 
electromagnetic spin interaction ([72] ) is 

£gm — 7, — RabcdS a S C . (73) 

8m 
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One can easily show that expressions ( |73|) and (|66| ) coincide (to within a 
factor e/m). It is sufficient to this end to take into account the relation 
S ab = e abcd S c u d , as well as the identity 

-^abcd ^ ^afe ^cd ^^fgh -ttabcd j 

which is valid for a sourceless gravitational field. 
By analogy with the magnetic moment 

2m 

it is natural to define the gravimagnetic moment 

_ t^tab Qcd 

2m 

The gravimagnetic ratio re, like the gyromagnetic ratio g in electrodynam- 
ics, may have in general any value. However, according to our semiclassical 
arguments, the value re = 1 in gravity is as preferable as g = 2 in electrody- 
namics. In any case, at g = 2 and re = 1 the spin equations of motion have 
the simplest form. 



On the other hand, it has been shown in JO]] that just this value of the 
gravimagnetic ratio, re = 1, follows from the wave equations in the Feynman 
gauge for the photon and graviton in an external gravitational field. The 
same value, re = 1, follows from the squared Rarita-Schwinger equation for 
s = 3/2 in a gravitational field [IT]. 



A second-order wave equation for an arbitrary spin in a gravitational 
background had been proposed long ago in PHI - For integer spins, its form 



corresponds also to re = 1. However, the value of re proposed in [[Hp for 
half-integer spins is quite different, it does not tend to unity even in the 
semiclassical limit s — > oo. Clearly, such a prescription does not look reason- 
able: in the semiclassical limit s — >• oo the parameter re should not change 
when the spin changes by 1/2. 

The situation with spin 1/2 is rather special [|lTJ. The properties of the 
spin matrices for s = 1/2, 

b =^(7 7 -7 7 ), 
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are such that the gravimagnetic interaction ([73| ) degenerates here into the 
mere scalar curvature 

" R, (74) 



16m 

without any consequences whatsoever for the motion of the spin. Obviously, 
the coefficient at this spin-independent stucture cannot be fixed by the above 
arguments. Formula (|74|) would predict for k = 1 the following wave equation 
for s = 1/2 

(-g^D^ - m 2 + ±R\i/> = (75) 

(the factor l/2m should be deleted from the expressions (|73|) and (|74]) when 
going over from a Lagrangian in our normalization to a wave equation). 
Meanwhile, the squared Dirac equation in a gravitational field is 



-g^D„D u -m 2 + ^R)r= 0. (76) 



The discrepancy between (|75|) and ([7(]) can be removed, for instance, by 
adding 

16m 

to the Lagrangian (^) for half-integer spins. This term in no way influences 
our semiclassical arguments, vanishing in the limit h — > 0. As to the recent 
proposal |42| to cure the discussed discrepancy by ascribing to the electron 



(which has no gravimagnetic interaction at all) the gravimagnetic ratio k = 2, 
we cannot see in it any real physical meaning. 



7 Multipoles of Black Holes 

But let us come back from elementary particles to macroscopic bodies. For 
a classical object the values of both parameters g and k depend in general 
on the various properties of the body. However, for black holes the situation 
is different. It has been shown in from an analysis of the Kerr- Newman 
solution that the gyromagnetic ratio of a charged rotating black hole is uni- 
versal (and equal to that of the electron!): g = 2. 

We will show that for the Kerr black hole the gravimagnetic ratio is k = 1. 
This value follows in fact from the analysis of the motion of spin of a black 
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hole in an external field done in [51 (though this statement was not explicitly 
formulated there). We will present here an independent and, in our opinion, 
simpler derivation of this important result. 

At great distance from a Kerr hole, the hole can be considered as a point 
source of a weak gravitational field. To linear approximation in the field of 
a hole at rest, the Lagrangian density corresponding to the interaction (|66|) 
can be written as 

L = 7— (5V) 2 h 00 5(r) . (77) 

The thus induced correction to the energy-momentum tensor has a single 
component: 

5T 00 = -^(^) 2 ^)- (78) 

In the gauge 

h^,u= 0, V= V- \^vK (79) 

the static Einstein equation for the corresponding correction hoo to the 00- 
component of the metric is 

A /ioo = 8irkToo . 

The correction itself is 

h 00 = k- (sV) 2 - . (80) 
m r 

Let us compare /loo with the corresponding contribution to the Kerr met- 
ric. In the Boyer-Lindquist coordinates this metric is 

ds z = (i _ Z£\ di * _ I± dr 2 _ Y,d6 2 - (r 2 + a 2 + T -^- sin 2 6)r 2 sin 2 9 

+ ^£1 sin 2 Q d(j)dt t ( gl ) 

where A = r 2 — r g r + a 2 , E = r 2 + a 2 cos 2 # , a = s/m. At r 9 = the 
metric fl8"T| ) describes a flat space in spheroidal coordinates [J5H]. Meanwhile, 
it is Cartesian coordinates which correspond in the flat space to the gauge 
( |79|) . The transition from the spheroidal coordinates to Cartesian ones is 
carried out with the required accuracy by the substitution 

_ _ a(ar) — fa 2 

r — > r + . 

2r 2 
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In the Cartesian coordinates the spin- dependent part of the 00-component 
of the metric 

Soo=l-^ + ^(3cos^-l) 



obviously coincides with hoo from formula (|80|) at k = 1. Somewhat more 
tedious consideration of the space components of the Kerr metric leads to 
the same result, k — 1. 

Let us note that the motion of the Kerr black hole in an external gravi- 
tational field is not described by the Papapetrou equation even if one leaves 
aside the problem of spin-orbit interaction linear in spin. The point is that 



this equation refers to the case k = [[14]]. 

It is proven in the same way that for a charged Kerr hole as well the 
gravimagnetic ratio k = 1. Moreover, it can be proven that the electric 
quadrupole moment of a charged Kerr hole also equals 

Q = -2^, 

the value, at which the interaction quadratic in spin decreases with energy. 



It has been shown [44[] that other, higher multipoles of a charged Kerr hole 
as well possess just those values which guarantee that the interaction of 
any order in spin (but of course, linear in an external field) asymptotically 
decreases with increasing energy. 
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